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Abstract 

We present an initial analysis of the properties of the all-sky image obtained by the Far-Infrared 
Surveyor (FIS) onboard the AKARI satellite, at 65 /im (N60), 90 fim (WIDE-S), 140 fim (WIDE-L),and 
160 /im (N160). Absolute flux calibration was determined by comparing the data with the COBE/DIRBE 
data sets, and the intensity range was as wide as from a few MJy sr“^ to >1 GJy sr“^. The uncertainties 
are considered to be the standard deviations with respect to the DIRBE data, and they are less than 
10% for intensities above 10, 3, 25, and 26 MJy sr“^ at the N60, WIDE-S, WIDE-L, and N160 bands, 
respectively. The characteristics of point sources in the image were also determined by stacking maps 
centred on photometric standard stars. The full width at half maxima of the point spread functions 
(PSEs) were 63", 78", and 88" at the N60, WIDE-S, and WIDE-L bands, respectively. The PSE at the 
N160 band was not obtained due to the sensitivity, but it is thought to be the same as that of the WIDE-L 
one. 
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1. Introduction 

The Japanese infrared astronomical satellite AKARI 
(Murakami et al. 2007) performed an all sky survey at six 
infrared bands centred at 9, 18, 65, 90, 140 and 160 jam. 
The four longer wavelength bands were observed with the 
Ear-Infrared Surveyor (EIS), one of the two focal plane in¬ 
struments onboard (Kawada et al. 2007). The two shorter 
EIS bands (N60 centred at 65 //m and WIDE-S centred at 
90 /im) were observed by the SW detector, and the other 
two bands (WIDE-L at 140 /im and N160 at 160 /im) by 
the LW detector. The all sky survey was started in 2006 
May and covered ^98 % of the whole sky during its cold 
phase (Nakagawa et al. 2007), continuing until helium ex¬ 
haustion in 2007 Aug. 

The All-sky Bright Source Catalogue (BSC, version 1), 
one of the primary goals of AKARI, was released in 2010 
Mar (Yamamura et al. 2010). Although an all-sky sur¬ 
vey image was not one of the primary goals of AKARI, it 
is scientifically valuable because extended sources are not 
covered by the BSC. There are two previous All-sky Ear- 
infrared (EIR) surveys that were obtained by the IRAS 
(Neugebauer et al. 1984) and COBE (Boggess et al. 1992) 
satellites. The IRAS survey was made with two EIR pho¬ 


tometric bands centred at 60 and 100 /im, with a spatial 
resolution of ^ 4'. The COBE mission expanded the wave¬ 
length coverage to 240 /im by using the DIRBE instrument 
(Hauser et al. 1998), but the spatial resolutions were only 
^0.°7. Recently, Miville-Deschenes & Lagache (2005) cre¬ 
ated a new data set of IRIS (Improved Reprocessing of the 
IRAS Survey) by combining the IRAS and COBE data. 
The IRIS is referred in many fields, but it only covers the 
wavelengths at or shorter than 100 jam. Since the thermal 
radiation from interstellar dust, whose typical tempera¬ 
ture is ^30 K, has its peak at 100-200 /im, the longer 
wavelength data are important. The present AKARI all¬ 
sky image includes data at 65, 90, 140, and 160 /im, with 
spatial resolution of ^1.'5. Thus we have created an all sky 
survey map from the AKARI data set (Doi et al. 2009) 
to provide new legacy data which has a wide potential 
application for use by the scientific community in many 
studies. In this paper, we report the calibration processes 
and performance of the EIS all sky maps. 

2. Data 

In the all sky survey mode, AKARI scanned the sky in 
an ecliptic orbit, with the telescope always pointing away 
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from the earth. At an altitude of ^700 km, AKARI made 
the orbit around the earth in ^100 minutes, scanning at a 
constant speed of ^3.'6 s“^. We refer to the all sky image 
obtained in this way as the AFASS (AKARI Far-infrared 
All Sky Survey) image. 

The observing strategy is explained in Kawada et al. 
(2007). The detector was read out continuously with a 
constant sampling rate for each array, and reset to dis¬ 
charge the photocurrent at intervals of about 2 seconds 
nominally, or 0.5 seconds toward bright sky. Close to 
the Galactic plane, where these long reset intervals were 
expected to lead to detector saturation, the detectors 
were reset for each sampling (correlated double sampling: 
CDS), with reset intervals corresponding to about 26 and 
45 ms for the SW and LW detectors, respectively. In this 
paper we refer to the “Normal mode” for the data which 
were taken with long reset intervals (i.e., 2 or 0.5 seconds) 
to distinguish it from the “CDS mode” for the short reset 
intervals. 

The data were pre-processed using the AKARI pipeline 
tool called “Green-Box”. This pipeline applied basic cal¬ 
ibrations such as corrections of non-linearity and sensi¬ 
tivity drifts of the detector, rejection of anomalous data 
due to charged particle hit (glitch), saturation, and other 
instrumental effects as well as dark-current subtraction. 
Following this initial processing, another “image” pipeline 
optimised for imaging large spatial scale structures was 
applied. At first, this image pipeline applied a transient 
response correction to the detector signal and subtraction 
of zodiacal emission. Since the detector had a slow re¬ 
sponse, the detector output drifted along the scan causing 
some stripe patterns in the image. To remove this artificial 
stripe pattern, we applied destriping. This destriping pro¬ 
cedure is as follows. Firstly, we applied a two-dimensional 
Fourier transformation of the image. Secondly, we masked 
the regions corresponding to the stripe patterns along the 
scan direction in wave number space. Next, we interpo¬ 
lated the intensities from the surrounding regions into the 
masked regions. Finally, we applied an inverse Fourier 
transformation to obtain the destriped image. After the 
image destriping, we performed recursive deglitching and 
flat-fielding to produce the final processed image. Details 
of the image pipeline are presented in Doi et al. (2009); Doi 
et al. (in prep.). The final image is made with a pixel scale 
of 15". 

3. Calibrating the FIS Images 
3.1. Calibration Method 

The FIS had a cold shutter whose temperature was 
^2.5 K and an internal calibration lamp, which had been 
well-calibrated in the laboratory. However, the character¬ 
istics changed in-flight because of differences of the oper¬ 
ation mode, surrounding environment, and modified sen¬ 
sitivities due to charged particle hits. Thus, we needed 
to re-calibrate the FIS data using in-flight data. In the 
FIS BSC, the flux densities of point sources were cali¬ 
brated by measuring the photometric standard stars, the 
asteroids, and the planets Uranus and Neptune with well- 
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known flux densities. This calibration was made for point 
sources and the data reduction processes were highly opti¬ 
mized for extracting point sources, for which the optimal 
method differed from that needed to reduce the all-sky 
image (which included diffuse emission). In addition to 
the normal all sky survey observations, AKARI also had 
capability to carry out pointed observations toward se¬ 
lected objects/regions with improved sensitivity. In these 
pointed observations, the FIS intensity scales were deter¬ 
mined on the basis of the COBE/DIRBE data by the self- 
flat-fielding method (Matsuura et al. 2007; Ikeda et al. 
2012). The calibration factors were obtained at those re¬ 
gions where the emission from the zodiacal light and the 
infrared cirrus dominated for the SW and LW detectors, 
respectively. However, this method was applied only to 
the Normal mode data and the brightness range of the se¬ 
lected regions is of the order of 1-10 MJy sr“^, while the 
brightness reaches as high as >1 GJy sr“^ at the Galactic 
plane in the all sky survey. Therefore, we needed to ex¬ 
tend this method to include brighter emission and also the 
CDS mode data to be able to calibrate the AEASS image. 

Eirst of all, we estimated the expected brightnesses in 
the four EIS bands. We used the DIRBE Zodi-Subtracted 
Mission Average (ZSMA) data set as an absolute bright¬ 
ness reference because it was one of the most reliable, well- 
calibrated all-sky data sets at far-infrared wavelengths 
to date. The uncertainties in the absolute brightness 
were 10-14 % for the 60-140 jim. bands (Hauser et al. 
1998). Although there was another data set of IRIS, it 
was also calibrated with the DIRBE, and it does not have 
>100 /im bands. The expected brightnesses were calcu¬ 
lated as follows: Eirst, we created continuous Spectral 
Energy Distributions (SEDs) by interpolating the DIRBE 
data with a power-low function, for each DIRBE pixel. 
Since the DIRBE data were calibrated under the flat spec¬ 
trum assumption, we made a colour correction using the 
applied power-law spectra. Second, we multiplied the 
interpolated SEDs by the EIS relative spectral response 
function. Third, we again made a colour correction to 
acquire the expected EIS values under the flat spectrum 
assumption. Eurthermore, since we selected ZSMA data, 
we calculated the brightnesses of the zodiacal emission 
in the EIS bands at the corresponding epoch using the 
Gorjian model (Gorjian et al. 2000). Here, we note that 
we only considered the “smooth cloud” component be¬ 
cause remaining small-scale structures, such as asteroidal 
dust bands and a mean motion resonance were not well un¬ 
derstood so far (Ootsubo et al. in prep.). The sum of the 
“DIRBE-expected” and zodiacal emission model bright¬ 
nesses should be compared to the observed EIS brightness. 

To calibrate the EIS data, we also prepared the AEASS 
image. The Green-Box tool converted the detector out¬ 
puts to intensities in units of Jy pixel” ^ with the conver¬ 
sion factors obtained in the laboratory. As noted above, 
they were not well calibrated for the in-flight environment. 
Although the image pipeline basically subtracted the zo¬ 
diacal emission, we skipped this process to create image 
at this time. Eurthermore, since the intensity of the zo¬ 
diacal light varies with time, we needed to restrict the 
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epoch. The all sky survey was the top priority during 
the first half year of AKARFs cold operations, with only 
a limited number of pointed observations made. During 
the remaining time, the all sky survey was completed at 
high priority, although the fraction of time devoted to 
pointed observations was increased. Thus, we chose data 
from the first half year (i.e., 2006 May to Oct), where 
there were few observational gaps due to pointed observa¬ 
tions. We convolved the FIS image with the beam pattern 
of the DIRBE at each DIRBE pixel to directly compare 
the EIS and DIRBE-expected brightnesses. Erom this 
comparison, we estimated the conversion factors from the 
Green-Box-processed data into the absolute intensity. As 
noted before, since we did not consider small-scale struc¬ 
tures of the zodiacal emission, we excluded those regions 
of 1/3 (Ecliptic latitude)! <40 deg to avoid the effects of 
zodiacal emission. The expected brightnesses of small- 
scale structures near the ecliptic plane are <5 MJy sr“^, 
<4 MJy sr“^, <1 MJy sr“^, and <1 MJy sr“^ for the N60, 
WIDE-S, WIDE-L, and N160 bands, respectively (Gorjian 
et al. 2000; Ootsubo et al. in prep.). After determining 
the conversion factors, we recreated the AEASS image us¬ 
ing the conversion factors, and then compared with those 
predicted from the DIRBE data again. This process was 
repeated several times, because the image pipeline does 
not guarantee the linear relation between the Green-Box 
calibrated EIS data and the processed image. Eigure 1 
shows the converged conversion factors. 

S.2. Flux Uncertainties 

Eigure 2 shows the pixel-to-pixel comparison between 
the calibrated EIS intensities and the expected values 
calculated from the DIRBE data, for each DIRBE pixel 
where |/3| >40 deg. In this figure, the zodiacal emission 
of the smooth cloud component was subtracted, and all- 
season data were used for the AEASS image. We note 
that the EIS data still included other zodiacal emission 
components of asteroidal dust bands and a mean mo¬ 
tion resonance. A good correlation exists for brighter 
(^10 MJy sr“^) region, while there are some deviations 
from the linear relation at fainter brightness levels due 
to the incompleteness of the zodiacal emission model and 
the DIRBE sensitivity of ^10 MJy sr“^. The intensity 
uncertainties were considered to be the standard devia¬ 
tions with respect to the DIRBE-expected intensities in 
Eigure 2 for various intensity ranges, as shown in Table 1. 
We concluded that the standard deviations are less than 
10 % for >10, 3, 25, and 26 MJy sr“^ at the N60, WIDE- 
S, WIDE-L, and N160 bands, respectively. We note that 
these estimates are based on the spatial resolution of 
DIRBE, i.e., about 0.7 deg. The WIDE-S band is the 
most sensitive one, and the two longer wavelength bands 
have large scatter at fainter (a few MJy sr“^) regions. 
This is because the performance of the LW detector is not 
as good as that of the SW one, which can be also seen in 
Eigure 1. Since the LW detector consisted of independent 
Ge:Ga elements, the sensitivity between the elements was 
largely scattered, while the SW detector was a monolithic 
Ge:Ga array and the sensitivity is quite uniform. It is hard 


Table 1. Uncertainties in the intensity calibration for each 
band. The uncertainties from the slope of unity in Figure 2 
are given for individual intensity bins. The 1 cr uncertain¬ 
ties is calculated for each DIRBE pixel where |/3| >40 deg. 



N60 

WIDE-S 

WIDE-L 

N160 

MJy sr ^ 

2.5 

* 

15.1% 

* 

* 

4.0 

* 

8.1% 

* 

* 

6.3 

20 . 1 % 

4.7% 

* 

* 

10 

9 . 9 % 

4.5% 

45 . 4 % 

56.1% 

16 

8 . 4 % 

4.2% 

14.3% 

16.6% 

25 

9 . 7 % 

4.3% 

10.1% 

11.2% 

40 

10 . 8 % 

4.5% 

8.3% 

8.8% 

63 

11 . 1 % 

5.3% 

8.1% 

7.9% 

100 

12 . 5 % 

5.7% 

8.9% 

8.9% 

160 

10 . 4 % 

7.3% 

9.6% 

9.7% 

250 

13.2% 

8.3% 

9.8% 

10.3% 

400 

11.3% 

8.1% 

8.8% 

8.9% 

630 

11.7% 

10.3% 

7.0% 

7.3% 

1000 

13.0% 

7.9% 

5.9% 

6.5% 


to estimate the sensitivity for the N60 band, because it is 
the one most seriously affected by the zodiacal emission 
of the 4 EIS bands. Although we excluded the regions 
of 1/31 <40 deg, where the zodiacal emission is dominant, 
there still remain the small-scale structures of the emis¬ 
sion at 1/31 >40 deg. The faint intensity (^1 MJy sr“^) 
regions are thought to be affected by the zodiacal emis¬ 
sion. Therefore, we should re-calibrate the image when 
the model of the zodiacal emission becomes more reliable. 

4. Point Sources 

4 . 1 . Point Spread Function 

Arimatsu et al. (2014) have investigated the character¬ 
istics of point sources in the EIS maps using preliminary 
data, which were almost the same as the released data. 
Here, we re-estimated the Point Spread Eunctions (PSEs) 
using the released data, following Arimatsu et al. (2014). 
We chose 352 photometric standard stars from Cohen et 
al. (1999) with expected flux densities of 0.02-10 Jy at 
the WIDE-S band. To derive the PSE, we picked up those 
sources whose flux densities were greater than the 2.5 a 
detection limit to secure sufficient number of sources. We 
then stacked 49, 97, and 9 stars with flux densities above 

1.2, 0.28, and 0.7 Jy at the N60, WIDE-S, and WIDE- 
L bands, respectively. The derived PSEs are shown in 
Eigure 3 and their full width at half maximum (EWHM) 
are listed in Table 2. We did not determine the PSE at the 
N160 band, because there were not enough bright stan¬ 
dard stars detected at the band. However, since the EIS 
PSEs are determined by the detector characteristics rather 
than optics and the PSE shapes at the WIDE-L and N160 
bands are similar in pointed observations (Shirahata et 
al. 2009), we assume that the N160 PSE is same as the 
WIDE-L one in the AEASS image. 
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Fig. 1. Relations between the Green-Box-precessed FIS intensity and the expected one based on the DIRBE data, including the 
zodiacal light emission. The colour scales are the number of the DIRBE pixels, where |/3| >40 degrees. The top row panels are for 
the Normal mode data, and the bottom rows for the CDS mode. The red broken lines show the adopted conversion factors. 
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Fig. 2. Relations between the calibrated FIS intensity and the expected intensity based on the DIRBE data, after subtraction 
of the smooth cloud component. The colour scales are the number of the DIRBE pixels, where |/3| >40 degrees. The broken lines 
indicate the deviations of ±3 MJy sr“^ from the linear relation of (expected intensity) = (calibrated intensity). 



Fig. 3. Obtained PSFs for the FIS N60 (65 fim), WIDE-S (90 /xm), and WIDE-L (140 i^m) bands. The contours are at 75 %, 
50 %, 25 %, and 10 % of the peak intensities. 











Calibration of AFASS data. 


5 


No. ] 


Table 2. FWHMs of stacked PSFs. 



N60 

WIDE-S 

WIDE-L 

EWHM [arcsec] 

63.4 ±0.2 

77.8 ±0.2 

88.3 ±0.9 

(in-scan) 

76.1 ±0.4 

102.3 ±0.3 

98.3 ±1.4 

(cross-scan) 

31.3±0.3 

55.0±0.1 

72.1 ±1.1 

solid angle [arcsec^] 

2699 ± 40 

6375 ± 30 

8031 ±237 

hlsource 

49 

97 

9 


^.2. Flux Calibrations 

We examined the flux calibrations for the point sources 
in the AFASS map. We prepared 10' x 10' maps centred 
at the positions of selected standard stars. We subtracted 
the sky background from each map by taking the median 
value, and normalized it to the expected flux density. We 
separated selected sources into several groups over their 
appropriate flux density. Then, we stacked the maps for 
each group and took average values. We performed aper¬ 
ture photometry on the stacked maps using an aperture 
of 90" radius with a sky background area of 120-300" 
radius. Figure 4 shows the observed-to-expected flux den¬ 
sity ratio for the stacked sources as a function of the ex¬ 
pected flux density in each bin, and the mean ratios are 
0.627T0.029, 0.696T0.008, and 0.381T0.043 for the N60, 
WIDE-S and WIDE-L bands, respectively. Due to sev¬ 
eral improvements of the data reduction processes, these 
ratios differe from those of Arimatsu et al. (2014), which 
used preliminary data. One significant difference is that 
the ratio at the N60 band becomes smaller for the faint 
source bins (<0.5 Jy). However, these small ratios are 
likely caused by the poor sensitivity of the N60 band, i.e., 
the signal-to-noise ratios of the faint sources in each map 
are smaller than one. Details of these processes are de¬ 
scribed in Arimatsu et al. (2014). 

4 . 3 . Comparison with the FIS Bright Source Catalogue 

We checked the confidence of the data by comparing 
with the EIS BSC. Here, we used sources whose qual¬ 
ity flags (EQUALxx) are 3 in the BSC, indicating the 
sources are reliable. We performed centroid determina¬ 
tion for each BSC source using the GCNTRD procedure 
in the Interactive Data Language (IDL) Astronomy User’s 
Library (Landsman 1993) to check whether or not the 
source is detected in the AEASS map. We rejected those 
sources whose positional differences between the BSC and 
the AEASS maps were larger than 48", which was used 
to extract a unique source in the BSC. As a result, we 
have detection rates of ^50 % for all the 4 bands. These 
relatively low detection rates were caused by the large 
size of the PSEs of the AEASS maps. Since most of the 
undetected sources are located around extended objects 
such as molecular clouds, planetary nebulae and galax¬ 
ies, it is hard to determine the source centroid in such a 
region (see Eigure 5). We note that the data reduction 
processes of the BSC were highly optimized for extracting 
point sources. Eigure 6 shows the positional accuracy of 
the AEASS maps. Eor the N60 and WIDE-S bands, the 



expected flux [Jy] 


Fig. 4. Observed-to-expected flux density ratio as a function 
of the expected flux density for the FIS (a) N60, (b) WIDE-S, 
and (c) WIDE-L bands. The error bars on the y-axis represent 
the 1 cr photometry uncertainties. The solid lines with the 
fllled areas are the weighted averages with 1 cr uncertainties. 
The dashed lines indicate the ratio of unity. The vertical lines 
represent the 5 cr detection limits of 2.4, 0.55, and 1.4 Jy for 
the All-Sky Survey at the FIS N60, WIDE-S, and WIDE-L 
bands, respectively (Kawada et al. 2007). 

accuracy becomes worse along the y-axis, i.e., the scan di¬ 
rection, caused by the elongated shape of the PSEs along 
the scan direction. 

We also performed aperture photometry for the de¬ 
tected sources using an aperture size of 90" with a sky 
region of 120-300" in radius. Since the aperture size was 
very large, we only considered the source observed in the 
Normal mode to avoid contaminations from surroundings. 
The CDS mode was used near the Galactic plane, and 
thus, the source density was high. Eigure 7 shows the ra¬ 
tio of the measured to BSC flux densities as a function 
of NDENS, which is the number of nearby point sources 
within 5' from each detected source in the BSC. The ratio 
is found to become large at NDENS>1, because of serious 
contamination caused by the large aperture size. 

We also show the flux density ratio as a function of 
the BSC flux density for the sources with NDENS=0 in 
Eigure 8. There is no significant correlation with the 
source flux density, which is consistent with Arimatsu et 
al. (2014). The large scatters in flux density ratio are due 
to the difference in photometry method, i.e., the BSC flux 
densities are determined by the PSE-fitting method. 
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Fig. 6. Positional differences of the point sources between the AFASS maps and the FIS BSC for the Normal (top rows) and CDS 
(bottom rows) modes. The dash-dotted ovals represent the FWHM of the PSF, tabulated in Table 2. The differences are estimated 
on the Ecliptic coordinates, and thus, the x- and y-axis correspond to the cross-scan and in-scan directions of AKARI, respectively. 


N60 WIDE-S WIDE-L N160 



NDENS 


Fig. 7. Flux density ratios between the AFASS map and the FIS BSC as a function of the NDENS number. The dotted lines 
show the ratio of unity, and the broken lines indicate the ratios derived from the PSF analysis (0.627, 0.696, and 0.381 for the N60, 
WIDE-S, and WIDE-L bands, respectively). 
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Fig. 8. Flux density ratios of the NDENS=0 sources between the AFASS map and the FIS BSC as a function of BSC ffux density. 
The dotted and broken lines have the same meaning as in Figure 7. 
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Fig. 5. Example of the centroid determination of BSC 
sources around M33 on the AFASS map in the WIDE-S band. 
The circle symbols are the sources whose centroids were de¬ 
termined, and the cross ones not well determined. 


Table 3. Number 
sources detected 

of 

in the 

the 

AFASS 

point 

maps. 

N60 

WIDE-S 

WIDE-L 

N160 

Normal mode 
BSC 10799 

AFASS 8683 

250804 

166831 

66276 

39415 

14722 

9353 

CDS modeT 

BSC 17980 

AFASS 12774 

122749 

50683 

52983 

23139 

22135 

11547 


(’!'Sources which were took in the CDS mode are indicated 
by the first bit of the FLAGSxx flag in the BSC.) 

5. Summary 

This paper presents an initial analysis of the imag¬ 
ing performance of the AKARI Far-infrared all sky sur¬ 
vey. We have performed the absolute calibration of the 
AFASS image by comparing with the COBE/DIRBE 
data, and the standard deviations with respect to the 
DIRBE-expected intensities are less than 10% for intensi¬ 
ties above 10, 3, 25, and 26 MJy sr“^ at the N60, WIDE-S, 
WIDE-L and N160 bands, respectively. The accuracy at 
the N60 band was reduced by the unreliability of the zodi¬ 
acal light model, and those at the two longer wavelengths 
were due to the detector performances. The AEASS im¬ 
age is improved in spatial resolution and wavelength cover¬ 
age, compared with the COBE/DIRBE and IRAS surveys. 
The positional accuracies are in good agreement with the 
PSEs. The characteristics of point sources in the image 
were also checked, and are consistent with Arimatsu et al. 


(2014). 

The AEASS data will be useful for studies of the Solar 
system, the Galaxy and nearby galaxies, in terms of both 
the spatial distribution of dust and their spectral en¬ 
ergy distributions. It is also helpful for future/on-going 
infrared/submm missions such as ALMA, JWST, and 
SPIGA. 

This research is based on observations with AKARI, a 
JAXA project with the participation of ESA. This work 
has been supported by JSPS KAKENHI Grant Number 
19204020, 21111005, 25247016, and 25400220. 
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